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Abstract 



We propose a new scenario for the bouncing universe in a simple five- dimensional 
braneworld model in the framework of Einstein-Gauss-Bonnet gravity, which 
works even with ordinary matter on the brane. In this scenario, the so-called 
branch singularity located at a finite physical radius in the bulk spacetime plays 
an essential role. We show that a three-brane moving in the bulk may reach 
and pass through it in spite of the fact that it is a curvature singularity. The 
bulk spacetime is extended beyond the branch singularity in the C° sense and 
then the branch singularity is identified as a massive thin shell. From the bulk 
point of view, this process is the collision of the three-brane with the shell of 
branch singularity. From the point of view on the brane, this process is a sud- 
den transition from the collapsing phase to the expanding phase of the universe. 
The present result opens a completely new possibility to achieve the bouncing 
brane universe as a higher-curvature effect. 
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1 Introduction 



Initial singularity is one of the most serious problems in modern cosmology. It is inevitable 
under several physically reasonable conditions in general relativity as a consequence of the 
singularity theorems [H [2] . Since the quantum effect of gravity dominates close to curva- 
ture singularities, the completion of the quantum theory of gravity is necessary to answer 
to this problem. Although the full quantum theory of gravity is not available at present, 
(Super)string/M-theory is one of the most promising candidates. Given the present circum- 
stances, the results in the low-energy effective theory must be useful to gain insights and 
suggestions for the problem. As a possibility to avoid the initial singularity, the bouncing 
universe is one intriguing scenario |3j. Because the effective energy- momentum tensor given 
from the semiclassical correction may violate the energy conditions, the quantum-gravity 
effect could cause the big bounce of the universe, a transition from the contracting phase 
to the expanding phase. (See [1] for a review.) 

Higher- dimensional cosmology based on string-generated gravity models is an ambitious 
way to explain problems in modern cosmology such as dark energy or initial singularity in 
a harmonic way by the effect of extra dimensions. After the publication of the papers by 
Randall and Sundrum [5], string- inspired braneworld cosmology, in which the observable 
universe is a (3+l)-dimensional timelike hypersurface (three-brane) embedded in a higher- 
dimensional anti-de Sitter (AdS) bulk space-time, has been investigated enthusiastically. 
(See [6j for a review.) The dynamics of the Friedmann- Robertson- Walker (FRW) brane 
universe has been fully investigated to explain cosmic inflation in the early universe or 
dark energy. In the present paper, we focus on the possibility of the bouncing brane 
universe. 

Among two Randall- Sundrum models, we consider the single-brane model. Assuming a 
perfect fluid obeying an equation of state p = (7 — l)p (0 < 7 < 2) on the FRW brane in 
the simplest braneworld with Z2 symmetry [7j, we can show that the bouncing universe is 
not realized for 2/3 < 7 < 2. This implies that, although the condition is milder than the 
standard four-dimensional FRW cosmology, a matter field with negative pressure violating 
the strong energy condition is still necessary for the big-bounce. Moreover, the bouncing 
universe requires the negative mass parameter in the Schwarzschild-Tangherlini-type bulk 
solution which causes a naked singularity. 

Nevertheless, there are models providing the bouncing universe on the brane with ordinary 
matter. One possible way is to introduce a matter field in the bulk spacetime. The 
low-energy effective theory of (Super)string/M-theory predicts several fields in the bulk 
spacetime such as a dilaton field or gauge fields. For example, by the dimensional reduction 
from the 11- dimensional Hofava-Witten model to the five-dimensional bulk spacetime, a 
dilaton field and a U(l) gauge field appear [S]. In the presence of a U(l) gauge field in 
the bulk, the bouncing universe is generically realized on the brane [9]. However, it is 
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expected that the bulk spacetime suffers from the mass inflation instability [TOl [TT] when 
the brane approaches the inner Cauchy horizon in the bulk [12] . In this research direction, 
a model of the bouncing brane universe without a bulk inner horizon was constructed in the 
presence of a SU(2) Yang-Mills field, in which the bulk spacetime is regular and free-from 
singularities [131 HI] ■ 

In the present paper, we reconsider the possibility of the bouncing brane universe with 
ordinary matter on the brane. As mentioned above, this is impossible in the standard 
Randall- Sundrum scenario based on general relativity unless the extra dimension is space- 
like. Actually the bouncing universe is generically realized if the extra dimension is time- 
like [15], but we do not consider this radical possibility here. Instead, we consider the 
braneworld in the framework of Einstein- Gauss-Bonnet gravity. While the bouncing uni- 
verse in the Gauss-Bonnet braneworld was discussed also in [TH] with a U(l) gauge field in 
the bulk, we consider the vacuum bulk here. 

Einstein-Gauss-Bonnet gravity is a natural higher- dimensional generalization of general rel- 
ativity as a quasilinear second-order theory. Also, the quadratic Lanczos (Gauss-Bonnet) 
Lagrangian appears in the low-energy limit of heterotic string theory together with a dila- 
ton [T7]. The early stage of the research history of the Gauss-Bonnet braneworld was 
somehow winding. There was a disagreement among the authors on the junction condition 
for the bulk spacetime and hence the resulting Friedmann equation on the brane. This 
confusion was settled finally by the establishment of the generalized Isreal junction con- 
dition by Davis [18] and independently by Gravanis and Willison [19]. (See also [20] and 
the comment in [H].) The correct Friedmann equation in the Gauss-Bonnet braneworld 
was first derived by Charmousis and Dufaux in [2T] and the dynamics of the Friedmann 
brane has been investigated by many authors up to now [221 [23] . (See [23] for the covariant 
gravitational equations on the brane.) Also, its viability [25], creation of the universe |26j . 
cosmological perturbations [27] have been investigated. The maximally symmetric brane 
was studied in p8]. Other aspects of the Gauss-Bonnet braneworld have also been stud- 
ied [21]. (See Sec. 5.7 in [30] for a review.) The purpose of the present paper is to point 
out a new and intriguing dynamical property of the Gauss-Bonnet brane universe. 

In Einstein- Gauss-Bonnet gravity, there is an exact symmetric vacuum solution which is 
a generalization of the Schwarzschild-Tangherlini solution in general relativity [311 132] . 
One of the most characteristic properties of this solution is the existence of a curvature 
singularity at a nonzero physical radius for negative mass parameter, which is called the 
branch singularity. The qualitative behavior of the Gauss-Bonnet brane universe has been 
analyzed by many authors; however, the role of the branch singularity in this context has 
not been clarified yet. In the general relativistic case, the moment when the brane reaches 
the central singularity in the bulk corresponds to the big-bang (or big-crunch) time on the 
brane. Therefore, one might naively think that the brane universe ends in (or starts from) 
some curvature singularity on the brane when it hits the branch singularity in the bulk. 
However, we will show that it is not the case. 
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In the present paper, we will show that the branch singularity is harmless in the sense 
that a finite body is able to reach there without being crushed to a point or ripped apart. 
We will also show that the vacuum bulk spacetime can be extended beyond the branch 
singularity in the C° sense and then the branch singularity is considered as a massive thin 
shell. As a consequence, it is concluded that the brane may reach and pass through the 
branch singularity into the extended region of the spacetime. This process provides a new 
scenario for the bouncing universe on the brane. 

The rest of the present paper is constituted as follows. In the following section, we present 
the bulk spacetime and the dynamical equation for the FRW brane in Einstein-Gauss- 
Bonnet gravity. In Sec. Ill, we explain our new scenario for the bouncing universe. In 
Sec. IV, the weakness of the branch singularity is shown. Concluding remarks and discus- 
sions including future prospects are summarized in Sec. V. In Appendix A, we present the 
geometric quantities of the bulk spacetime. In Appendix B, the asymptotic behavior of the 
brane universe for a — )■ cxd is presented. 

Our basic notation follows [33]. The convention for the Riemann curvature tensor is 
[Vp, Vo-]V^ = R^upaV^ and Rfj,u = R^^pu- The Minkowski metric is taken as diag(— , +, +, +, +), 
and Greek indices run over all spacetime indices. We adopt the units in which only the 
five-dimensional gravitational constant G5 is retained. 



2 Preliminaries 

We consider the following five- dimensional action for the bulk spacetime: 



where A is the cosmological constant and K5 is defined by the five-dimensional gravitational 
constant as := y^SirG^. The Gauss-Bonnet term Lgb is a combination of the Ricci 
scalar R, the Ricci tensor and the Riemann tensor R^upa as 



which does not give any higher-derivative (more than the second-derivative) term in the 
field equations. The constant a in the action is the coupling constant of the Gauss-Bonnet 
term and for a ^ our model reduces to the Randall- Sundrum model. 

The gravitational equation given from the action (12. ip is 





■— R^ — AR^yR^" + R^ypaR'^"'"^ ^ 



(2.2) 



V 



(2.3) 
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where 



G^u ■=R^iu — -^dfiuR, (2.4) 

Hfj,i, ■.=2{^R^y — 2R^aR°'v — ^R"^ Rfj,aul3 + Rfi"'^^ Rual3-f^ — -Q^iuLqB- (2.5) 

The Gauss-Bonnet term in the action is obtained in the low-energy hmit of heterotic su- 
perstring theory together with a dilaton in ten dimensions [T7], in which case a is regarded 
as the inverse string tension and positive definite. We therefore assume a > throughout 
this paper. We also assume A < and 1 + 4q;A/3 > in addition, the latter of which 
ensures the existence of nondegenerate maximally symmetric vacuum solutions. 



2.1 Bulk solution 



In this system, a vacuum solution is obtained as a warped product manifold M.^ x , 

where is a three-dimensional space of constant curvature. In the equations which follow, 
k denotes the curvature of and takes the values 1 (positive curvature), (zero curvature), 
and —1 (negative curvature). The metric of the vacuum solution is given by 



d/S ^ — Q 



- Hr)dt' + J^^+r' [dx' + fk{xf{d9^ + sin^ Odcf")] , (2.6) 



Kr):=k+'^{lT\jl + ^ + \aK\, (2.7) 



where /x is constant, /o(x) •= Xi fiix) '■= sinx, and /-i(x) := sinhx [311 [32]. The relation 
between /z and the global mass parameter M is 

= 7^, 2.8 

where the constant is a unit volume of K^^ if it is compact. In the asymptotically fiat 
case {k = 1), M gives the ADM mass. 

It is seen that the solution has two branches corresponding to the sign in the metric function 
h{r). We call the family with the minus (plus) sign the GR branch (non-GR branch). Only 
the GR branch solution has the general relativistic limit as 



lim h(r) = k - - -Ar^. (2.9) 

a-s-o 8r^ 6 



The maximally symmetric vacuum in the non-GR branch was shown to be unstable 
and so we only consider the solution in the GR branch in the present paper. 
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The global structures of this spacetime depending on the parameters have been clarified |35j . 
In this spacetime, there are two classes of curvature singularity for /i 7^ 0. One is the central 
singularity at r = and the other is the branch singularity at r = rb(> 0), where the term 
inside the square-root in the metric function (12.71) vanishes, rb is explicitly given by 

.,:=(-^V''. (2.10) 



3 + 4aAy 

The branch singularity exists if fi is negative. The metric and its inverse are finite at r = rt 
(but their derivatives blow up) and the metric becomes complex and hence unphysical at 
r < rb. 



The relation between the mass parameter n and the radius of the Killing horizon [defined 
by /i(rh) = 0] is given by 

fx = 16ae + 8krl - ^Ar^ =: /ih(rh). (2.11) 
The relation between fi and rb is given by 

fi = ~(^l + ^c^^yt =■■ /"b(rb). (2.12) 

The number of horizons, their properties, and the existence of singularities depending on 
the parameters are understood from the functional forms of fx = fxi^{r) and fx = fXh{r), 
which are drawn in Fig. [1] (See [35] for the detailed analysis.) We summarize the global 
structures as follows. For k = 1, there is one nondegenerate outer horizon for fi > fio := 
16a and no horizon for fi < fio- For k = 0, there is one nondegenerate outer horizon 
for fi > and no horizon for fx < 0. For k = —1, there is one nondegenerate outer 
horizon for fi > fXc{'-= — 16q;(1 + 4aA/3)), two nondegenerate (outer and inner) horizons 
for /icx(:= 12(1 +4aA/3)/A) < fx < fic, one degenerate horizon for fi = fx^x, and no horizon 

for /i < /iex- 



Near r = 0, the metric function in the GR branch behaves as 

which is valid for fx > 0. Therefore, the central singularity is timelike, null, and spacelike 
for k - {\l^\[Jifa > 0, k - {l/A)^J fi/a = 0, and k - {Ij^^fjxfa < 0, respectively. On 
the other hand, near the branch singularity, the metric function behaves as 

h{r) ^{k + ^^ -^^l + ^«A(r - r^Y^'. (2.14) 

Therefore, the branch singularity is timelike and spacelike for k + r^/ (4a) > and k + 
r^/(4a) < 0, respectively. The Penrose diagrams in the case of < 0, which is of our 
interest, are drawn in Fig. |21 
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|X k=l k=0 k=-l 




Figure 1: The curves fi = fih{r) and /i = fih{r) in the GR branch with A < 0, a > 0, and 
1 + 4a A/3 > 0. Thin curves correspond to /i = /ih(r) for each k. A thick curve corresponds 
to /i = fih{f)- The metric becomes complex and unphysical in the region of /i < fihif)- The 
curve /i = fih{r) with k = —1 terminates on the curve ^ = fih{r) at /i = /ic- 



2.2 Friedmann equation on the brane 



We consider a three-brane in the bulk spacetime fl2.6p . which is a timlike hypersurface 
described by r = a(r) and t = T{t), where the parameter r is the proper time on the 
brane. The tangent vector to the brane is written as 

T^ + d^, (2.15) 



dxf^ dt dr ' 

where a dot denotes the differentiation with respect to r. The normalization condition 
u^u^ = — 1 leads to 

and the induced metric of the three-brane cjab is given by 

dsl = Qabdy^dy^ = -dr^ + a(r)2 [dx^ + fk{xf{de^ + sin^ Odcj)^)] . (2.17) 
This is the FRW metric with the spatial curvature k. 

The dynamics of the three-brane, namely the behavior of the scale factor a(r) on the 
brane, is determined by the junction condition. Here we simply assume the Z2-symmetry 
of reflection with respect to the brane; we take two copies of the bulk spacetime with 
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Figure 2: The Penrose diagrams of the bulk spacetime (12.61) in the GR branch with a > 
0, A < 0, 1 + 4a; A/3 > 0, and yU < 0. A zig-zag hne and a dashed hne represent a 
branch singularity and infinity (3), respectively. A thin line and a double line represent a 
nondegenerate and a degenerate Killing horizon, respectively. The diagram (a) corresponds 
to the cases of A; = 1 and 0. For k = —1, the diagrams (b), (c), (d), and (a) correspond to 
the cases of /Xc < < 0, /Zex < < /^c, = f^ex, and fj, < fj,ex, respectively. 
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r < a(r) and paste them at r = a(r). The junction condition in Einstein- Gauss-Bonnet 
gravity is given by 



[K\]^ - 5\[K]^ + 2a(?,e[J\]^ - e5\[J]^ - 2P\AK''%) = -e^S^, (2.18 



where 

J,, : = i {2KKadK\ + K,fK^fK,, - 2K,,K^fKf, - K'K,,) , (2.19) 

Padbf ■=Radbf + 2ha[fRb]d + 2hd[bRf]a + Rha[bhf]d- (2.20) 

Parfb/ is the divergence-free part of the Riemann tensor, i.e., 

D,P%^,f = 0, (2.21) 

where Da is the covariant derivative on the brane. We have introduced the notation 

[X]±:=X+-X-, (2.22) 

where is the quantity X evaluated either on the + or — side of the brane. e = 1 and 
e = —1 are used for timehke and spacehke branes, respectively. For the dynamics of our 
three-brane, we take e = 1. We consistently assume the form of the energy-momentum 
tensor S'^h on the brane as 

5"% = diag(-p,p,p,p) + diag(-cr, -a, -a, -a, ), (2.23) 

where p and p are the energy density and pressure of a perfect fluid on the brane, and the 
constant a is the brane tension. 

The junction condition f l2.18p with e = 1 gives the modified Friedmann equation for the 
brane universe as 



3 V CL 



2 



(2.24) 



where H := a/ a |T8]. The energy-conservation equation on the brane is obtained in the 
standard form as 

p=-3H{p + p). (2.25) 

In order to close the system, we have to introduce an equation of state for matter. We 
assume the following linear equation of state: 

P = {l- 1)P- (2.26) 
Equation (12.25^ is then integrated to give 

P = % (2.27) 
9 



where a constant po is assumed to be positive, so that p is a monotonically decreasing 
function of a for 7 > 0. Because 7 = is equivalent to the cosmological constant, we do 
not consider this case. We assume < 7 < 2 which satisfies the dominant energy condition. 

In the Randall-Sundrum model, the Friedmann equation on the brane [7] becomes 

H' = i(^ + -X-^ + L + l^=-- VoM. (2.28) 



36 \a^^ J a2 8a^ 6 

Now we have a one- dimensional potential problem with one dynamical degree of freedom 
a(r) and the qualitative behavior of a is completely understood by the form of the potential 
Vgr. The region with VGR(a) > is the allowed region for dynamics. The above equation 
fl2.28p is compared with the standard four-dimensional Friedmann equation: 

H' = ^^--, + \Aa =■■ V;tand(a), (2.29) 

where is defined by the four- dimensional gravitational constant as /€4 := \/WkGI and 
A4 is the four- dimensional cosmological constant. 

In the presence of the Gauss-Bonnet term, the dynamics of the brane is drastically modified. 
Equation f l2.24p in the GR branch is written as 

^''»+,V=fA+i^+tfA(A+i±^+,^A\ (2.30) 



256^2 \^a37 J \a? 4a \a? 8a 

A: = l + ^ + ^aA. (2.31) 

Qj O 

We note that A < 1 is satisfied if /x < 0. The necessary conditions for the above equation 
to be physical are 



k 1 - VA , 
^>0, —+ / +H^>0. 2.32 

For our purpose, we rewrite Eq. f l2.30p in the form of if^ = VGB(ct)- Since Eq. fl2.30p is 
a cubic algebraic equation for H^, there can be three real roots at most. However, it is 
shown that there is only one positive real root at most. 

The proof is given as follows. For this purpose, we define a cubic function F{x) as 



2 



k , 1-y/A Wsfca 2 + y/A 

la ^^J I 3a2"^ 3 ' 3^"^/ 36 V 



4 / \ 2 

4 / Po ^ 



P(^)--=\-2+—^ + ^ \tI^ + ^^+o^^\ + (2-33) 



The Friedmann equation fl2.30p is now F{x) = with x = H^. The dynamics is allowed 
only if F{x) = admits real and positive roots. We show that F{x) = has only one real 
positive root at most. 
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We obtain 



dF _ 64a2 / ^ 2 + 
dx 3 " " 



k 2 
+ - 



A 



+ X 



and the solutions of the algebraic equation dF/dx = are x = x± (x+ > where 

-2±^/A k 



x± 



8a 



i2- 



Since we obtain 



F{x, 
Fix^ 



^3/2 



18« 36 V a?^ 



/ Po 



/ Po 



36 V 



a 



(2.34) 

(2.35) 

(2.36) 
(2.37) 



which satisfy F{x^) < F{x^) < 0, it is concluded that there is one real positive root for 
F{x) = at most. A sufficient condition for the existence of a real positive root is -F(O) < 0. 

For a > 0, Eq. ( ]2.30p is solved to give the following modified Friedmann equation on the 
brane: 



H'=Vc 



GB(+)l 



Vr 



GB(+)l 



1 



Ska 



2 + I _^ 256a^p2 ^ 16 J2a3p2 (^uga'^P^ + A^/^ 



(2.38) 

1/3 



+ a| ^3/2 + 256a^P^ + 16 W2a3p2 (^uSa^P^ + A^/^^ | 



-1/3- 



P2: = 



256^2 \ q37 



(2.39) 



(2.40) 



Here VGB(+)(a) denotes the effective potential for a > 0. It is clear that inside both cubic 
and square roots in the potential are non-negative. The dynamics is allowed in the domain 
of a where Vgb(+)('2) is real and positive and the conditions (12.321) are satisfied. Indeed, the 
second necessary condition in Eq. (12.321) is always satisfied. Using the form of the potential 
VGB(+)(a), we can write it in the following form: 

k 1-^/A 



=i_ 1^3/2 ^ 256a^P^ + 16^2a3p2 (^uSa^P^ + ^3/2^ | 



-1/3 



X 



-VA + 1^3/2 ^ 256a^P^ + 16W 2a3p2 (^uSa^P^ + A^/^^ | 



1/3-1 2 



=: WM 



(2.41) 
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It is obvious that this condition holds only for a > 0. 

The effective potential Vgb(-)(«) for a < is different from Vgb(+)('3') and given by 

1/3 



a 



Ska 



+ A\ -A^'^ - 25Qa^P^ + 16. / 2a^P^ f 128a3p2 + ^3/2 



(2.42) 



With this potential, the second necessary condition in Eq. fl2.32p becomes 



k l-VA 



^j-A^^/^ - 256a^p2 + 16W2a3p2|^i28a3p2 + ^3/2^ | 



-1/3 



+ - 256a^P^ + 164/2a3p2 |^128a3p2 + ^3/2^ | 



1/3- 



2 



(2.43) 



Therefore, the second necessary condition holds for a < 0. 

The condition that inside the square root in the potential is nonnegative is given by 

128a^P^ + A'/' = ^(P^ + a]\(l + ^ + -aA] "\ 0. (2.44) 
2 ya^'i' J \ ?, J 

This may be violated for some a. On the other hand, the condition that inside the cubic 
root in the potential is non-negative is given by 



A^'"^ + 256a^p2 <16y 2a3p2 |^i28«3p2 + ^3/2^ . (2.45) 

This is always satisfied under the condition ( I2.44p because the left-hand side is nonpositive. 
Since the equality in Eq. (I2.45p cannot be satisfied, inside the cubic root in the potential 
is positive definite. 

Let us see what happens if inside the square root in the potential becomes zero. There are 
two possibilities for that, P = and 128q;^P^ + ^3/2 = Q. It is shown that the derivative 
of the potential blows up only in the latter case. Therefore, the equality in Eq. (12.441) 
corresponds to a curvature singularity where a and a are finite but a blows up. 
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3 Bouncing Gauss-Bonnet braneworld 



The bouncing solution is characterized by the transition from the contracting phase (a < 0) 
to the expanding one (a > 0) of the universe. As seen in the previous section, once we 
write the Friedmann equation in the form of = V{a), V{a) > is the allowed region for 
dynamics. We assume that there is at least one domain of positive a with V{a) > since 
there is no dynamical solution otherwise. Then, the bounce (recollapse) in a conventional 
sense occurs at the lower (upper) bound of this domain a = qb satisfying V^qb) = 0. The 
evolution of the contracting universe in the domain a > ae momentarily stops at a = ob 
and then starts to expand in the domain a > ae- Therefore, the sufficient condition for 
this conventional bounce is ^(0) < 0. (Similarly, the sufficient condition for recollapse is 
lim^^oo V{a) < 0.) 

In this section, we discuss the possibility of the bounce by the asymptotic analysis for 
a — )■ 0. The asymptotic behavior for a — )■ oo is presented in Appendix B. We first consider 
the possibility of the conventional bounce in the Randall-Sundrum braneworld and the 
Gauss-Bonnet braneworld with /i > 0. We will see that, in comparison with them, the 
situation in the Gauss-Bonnet braneworld with < is very different. The primary reason 
is that a < Tb is not in the physical domain of a. We will see what happens when the brane 
hits the branch singularity in the bulk. 



3.1 Conventional bounce condition 



In this subsection, we consider the condition for the bounce in a conventional sense. In the 
standard Friedmann cosmology (12. 29 p . we obtain 

limK,..d(a)-|4^-4- (3-1) 

Hence, the bounce occurs only for /c = 1 where 0<7<2/3 or 7 = 2/3 with po < 3/^4 
is satisfied. In the case of /c = 1, 7 = 2/3 and po = S/kI, the spacetime is just (A)dS for 
A4(<) > 0. 

In the Randall-Sundrum braneworld (I2.28p . where we assume A < and /i 7^ 0, we obtain 

limyGR(a)^-4- + A- (3.2) 

Hence, the bounce occurs for < 7 < 2/3 with /i < 0. For 7 = 2/3, the bounce occurs for 
< /^GR{7=2/3), where 

/^GR(7=2/3) := (3.3) 
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If a precise relation between po and /i is satisfied as /x = /Ugr(7=2/3), the bounce condition 
is more complicated for 7 = 2/3. The bounce conditions are summarized in Tabled! 

Table 1: Conditions for the conventional bounce in the Randall- Sundrum braneworld for 
arbitrary k. 





< 7 < 2/3 


7 = 2/3 


2/3 < 7 < 2 


/i > 


No bounce 


No bounce 


No bounce 


fi <0 


Bounce 


(No) Bounce for fx < (>)/Ugr(7=2/3) 


No bounce 



In the Gauss-Bonnet braneworld f l2.38p . the bounce condition is modified a lot. Let us 
consider the case with /i > 0. Our purpose is to clarify the asymptotic behavior of the 
potential for a — )■ 0. For this purpose, the asymptotic form of W+{a) defined by Eq. fl2.4ip 
is useful. We write W+{a) in the following form: 



6aa 



{any/- 



4 . 
1 + -aA 



+ 256w^ + 



256 



VT28 



zW 



VTTl28uP 



-1/3 



1 + 256?^;^ + 



256 



VT28 



zW 



Vl + 128u;2 



1/3- 



=: L(w) 



w :=(«3^-3/2p2)i/2_ 
w is written in terms of a as 



w 



256(a/i)3/2 
which is expanded around a = as 



-aA 1 — 



-3/2 



1 + 



37 



Po 



1/2 



2„,,^4^|l/2^3-37 



lim w ~ , 



4 . \ 



era' 



37 



Po 



The derivative of w converges (diverges) forO<7<2/3 (2/3<7<2) near a = as 

a™ da ~ 16(a/i)3/4 
Hence, L{w) can be expanded around a = OforO<7<2/3as 



r r/ ^ ^12 2 
hm L UM ^ w 

a^o ^ ^ 9 



2pgaft:t 6-67 
9(«/i)3/2" ' 



from which we obtain 



lim W+(a) 

a->0 



,_Po^^4-67 



(3.4) 
(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 
(3.10) 



14 



Thus, asymptotic behavior of the potential around a = for 0<7<2/3is 



hm Vgb{+) ( 



a->0 



Aa J Aa 36q;/x 



(3.11) 



For 2/3 < 7 < 1, we only know lima_^o = 0{a^) and linia^o L{a) = 0{a^), where 6 and S 
are positive constants. Using them, we obtain 



liml^+(a) =0(a^-2), 

a— >-0 

limW-) = 4(*-^)~ + 0(. 



5-2\ 



(3.12) 
(3.13) 



from which we realize the leading term around a = 0. It is concluded from the analysis 
that, for < 7 < 1, no bounce occurs for /c = 0, — 1 or = 1 with fi > 16a(= /io), while 
the bounce occurs for A; = 1 with fi < fiQ. In the case of = 1 with = fiQ, the bounce 
occurs for < 7 < 2/3, while a more careful analysis is required for 2/3 < 7 < 1. For 
7 = 2/3, the (no) bounce condition is given as yU > {<)pIk'^/9. 



For 1 < 7 < 2, w blows up near a = and 

lini VGB(+)(a) 



Therefore, the bounce does not occur. 



(3.14) 



For 7 = 1, w behaves near a = as 



Using this, we obtain 



limV"GB(+)(a) 



X 



lim w — , 

a^O 16(a;U)3/4 



1 + 



Po 



4« 



+ 256wg + 



256 



1 + 256u;^ + 



256 
yT28~ 



7128 

l/3n 2 



(3.15) 



wo\/l + 128^2 
1 



-1/3 



4a 



+ 0(a). 



16(a;u)3/4 ■ 



(3.16) 



Hence, the bounce does not occur for /c = 0, — 1 because the coefficient of is positive. 
On the other hand, the bounce occurs (does not occur) for = 1 if the coefficient of is 
nonpositive (positive). The result is summarized in Table |2J 
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Table 2: Conditions for the conventional bounce in the Gauss-Bonnet braneworld with 
/i > 0. 







k = 1 


A; = 


k = -l 


< 7 < 2/3 


(No) 


bounce for /i < (>)/io 


No bounce 


No bounce 


2/3 < 7 < 1 


(No) 


bounce for /i < (>)/Uo 


No bounce 


No bounce 


7 = 1 




See the maintext 


No bounce 


No bounce 


1 < 7 < 2 




No bounce 


No bounce 


No bounce 



3.2 Novel bouncing Gauss-Bonnet braneworld for /i < 

We have seen the bounce conditions for /i > in the Gauss-Bonnet braneworld. In this 
subsection, we discuss the case with /x < 0, where the situation is drastically changed. As 
explained, a < rb is not in the physical domain of a and so we focus on the behavior of 
the potential around a = rb. If the potential is non- negative, the brane reaches a = rb. If 
the potential is negative near a = rb, the bounce (in the conventional sense) occurs or a 
singularity appears at some a > rb. 

First we present the condition that the brane hits the branch singularity in the bulk space- 
time. The function A defined by Eq. fl2.3ip behaves near a = rb as 

^(a)^-^(a-rb). (3.17) 

The behavior of Vgb(+)('^) around a = rb is given by 

V^GB(+)(a) ~rGB(+)(rb) + 0{a - rb), (3.18) 

where 

W)K)=-^-^ + P(rb)^/^ (3.19) 

Therefore, the brane reaches the branch singularity if VGB(+)(^b) > 0. This condition is 
satisfied if 

2 

h(r^) = k + ^<0. (3.20) 
4a 

Therefore, if the branch singularity is spacelike, which is realized only for k = —1, the 
brane reaches there. If the branch singularity is timelike, the condition VGB(+)('"b) > is 
written as 



37" 
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It is noted that this condition is not so sensitive about 7 and reahzed even with positive 
pressure (7 > 1). On the other hand, if 



is satisfied, a = rt is not in the physical domains of a and hence the bounce occurs or a 
singularity appears at some a > rb instead. 

We saw the condition under which the brane hits the branch singularity in the bulk. Let 
us see what happens then. Now the important fact is that the curvature invariants on the 
brane do not blow up even when the brane approaches the branch singularity in the bulk. 
In fact, the behavior of the scale factor a(r) near the branch singularity is obtained as 



where Tb is the cosmological time on the brane to reach the branch singularity. This is the 
Taylor series around t — t^, and hence the curvature invariants are all finite around there. 
However, because the allowed domain of a is a > rb, we must take the minus and plus signs 
of ai for r < Tb and r > Tb, respectively. As a result, the evolution near r = Tb represents 
the transition from the collapsing phase (r < Tb) to the expanding phase (r > Tb). Thus, 
the bouncing universe is realized on the brane. 

Prom the bulk point of view, this process is that the brane reaches the branch singularity 
and passes across it. It is emphasized that the spacetime on the brane is not free of 
singularities then, but there appears just a shell-type instantaneous singularity. In the 
generic case, the derivative of a(r) (velocity) is not continuous and the metric on the brane 
is C° at T = Tb. The junction condition on the brane then shows that there is a matter 
distribution on the spacelike hypersurface t = Tb on the brane. This means that a shell-type 
singularity appears instantaneously on the brane at t = Tb but it is rather harmless since 
it stems from the thin-shell approximation of the brane as well as the branch singularity. 
(In the next subsection, we will show that the branch singularity may be considered as a 
massive thin-shell.) With a fine-tuning giving ai = 0, in contrast, the metric on the brane 
becomes analytic around t = Tb and there is no shell-type singularity. 

Our claim is that the brane reaches the branch singularity and passes through it in the 
generic case and this process is the collision of the three-brane with the shell of branch 
singularity from the bulk point of view. In order to support this claim, we show in the next 
section that the branch singularity is indeed harmless for a finite body moving radially. 
Then, a next natural arises; what is the bulk spacetime on the other side of the branch 
singularity? In order to answer to this question, we have to consider the extension of the 
bulk spacetime beyond r = rb. 




(3.22) 




rb + ai(T-Tb) + 0((T-Tb)'), 
rb^GB(+)(rb), 



(3.23) 
(3.24) 
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3.3 C extension of the bulk beyond the branch singularity 

In this subsection, we discuss the extension of the bulk spacetime beyond the branch 
singularity. It is seen in fl2.6p that the metric and its inverse are finite at the branch 
singularity. This implies that the spacetime can be extended beyond r = rb in the C° sense. 
The extension is not unique in general; however, thanks to the Birkhoff's theorem in 
the system there are only two candidates for the extended spacetime, namely the GR 
and non-GR branches of the vacuum solution (12.61) . Among these two, the GR branch 
should be chosen because of its dynamical stability. In Fig. 3(a), the global structure of 
a C°-extended spacetime is drawn as an example. [It is the extension of the spacetime of 
Fig. 2(a) and the diagrams become different for the spacetimes of Figs. 2(b), (c), and (d).] 



(a) 



(b) 



Figure 3: The Penrose diagrams of (a) a C°-extended bulk spacetime beyond the branch 
singularity [corresponding to Fig. 2(a)] and (b) its regularized version with a dust fluid 
for k = —1 and /i < /igx- In (b), a dust fluid fllls in the shadowed region and the branch 
singularity is regularized and replaced by a regular wormhole throat (a thick solid line). 



We have constructed the C°-extended bulk spacetime. Now let us study the branch singu- 
larity in this spacetime in more detail. The induced metric on r = rt is given by 

dsl =gabdy''dy^ 

= - h{r^)de + rl [dx^ + MxYide' + sin' 9d^')] . (3.25) 

We know that the derivatives of h are not flnite at r = rb. The component of the extrinsic 
curvature Kat of the hypersurface r =constant diverges in the limit r — )■ rb. Nevertheless, 
the junction condition provides a flnite value of the energy-momentum tensor on r = rb. 
This implies that the branch singularity can be identifled as a massive shell in the extended 
bulk spacetime. (It is noted that, in Einstein-Gauss-Bonnet gravity, even vacuum 
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spacetimes can be constructed [36].) We consider the cases where r = r^, is timehke 
(/i(rb) > 0) and spacehke {h{r^) < 0), separately. (A more careful treatment is required in 
the case of h^r^,) = 0.) 

First we consider the case where r = rb is timelike. We have 

2 

h{rh) = k + 0) (3.26) 

and the tangent vector to a timelike hypersurface r =constant is written as 

n d 1 d , , 

which satisfies m^m^ = — 1. The unit normal one- form to this hypersurface is given by 

n^dx'' = j dr, (3.28) 

where n^u'^ = and n^n'^ = 1 are satisfied. 

The extrinsic curvature of this hypersurface is obtained from Kat, '■= (Vjyn^)e^e^, where 
:= We have 

e'jy'' = dt, eldy^ = 0, e^dy'^ = 5]dy\ (3.29) 

and the nonzero components of K'^^ are 

K\ = -^h'{r), K^^ = --^)5^,, (3.30) 
where a prime denotes differentiation with respect to r. 

Taking the limit r — )■ rt, A'^ blows up because h' blows up there. However, the junction 
condition fl2.18p (with s = 1) shows that S''^f, on the hypersurface remains finite for r — )■ rb 

as 

lim 5"% = diag(-pb,Pb,Pb,Pb), (3.31) 

where pb and pb are interpreted as the energy density and pressure at r = rb. Using the 
following fact, 

\iui\h'Uah{r)-Aak-r'^)\^ — ^, (3.32) 

r^rt J 2rg 



we actually obtain 



/i(rb)2/2 8/i(rb)r^-/x 
^ = -16a.— . = (3.33) 
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Thus, the branch singularity can be considered as a massive thin shell. Since pb < 0, the 
matter on r = rt violates the weak energy condition. 



Next we consider the case where r = rb is spacelike, which is realized only for /c = — 1. In 
this case, we consider the tangent vector to the (spacelike) hypersurface as 

' ' (3.34) 



which satisfies Uf^u'^ = 1. The unit normal one-form to the (spacelike) hypersurface is 
then given by 

n^dx^ = ^ dr, (3.35) 

-h(r) 



where ?t,„m^ = and = — 1 are satisfied. The nonzero components of K°-iy are 



K 



h'(r)_ K'j = -^/^h{^6'j. (3.36) 



Here we write the energy-momentum tensor S'^f, on the (spacelike) hypersurface as 

Jim S^b = diag(Pb(r), ^b(t), ^b(t), ^b(t)), (3.37) 

where Pb(r) and Pb(t) are radial pressure and tangential pressure at r = rb, respectively. 
Using Eq. f l3.32p . the junction condition fl2.18p with e = —1 gives 

(-/i(rb))3/2 8h(rM + f^ 

Pb« = 16a^^-J^, PHt) = A==^- (3-38) 

They are finite and hence the branch singularity may be considered as a massive spacelike 
thin shell. Since this is a spacelike shell, it is difficult to discuss the energy condition for 
the matter field there. 



3.4 A simple regularized bulk model with matter 

We have seen that the branch singularity is weak and a three-brane may pass through 
it. However, of course, the theory is no more valid around there in the bulk because it 
is a curvature singularity. To make matters worse, it is a naked singularity in most cases 
and so it causes a problem of the boundary condition when we consider, for example, the 
development of the perturbations in the bulk. However, the weakness of the singular- 
ity could suggest that it is regularized in the effective theory of the final theory. If so, 
the spacetime around the branch singularity should be described by some totally regular 
spacetime which approaches the vacuum spacetime in the far region. Here the correction to 
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the original theory around the singularity would appear in the field equation as an effective 
energy-momentum tensor, which might violate the energy conditions. 



In this subsection, we present an exact bulk model with a simple matter field as an example 
of such spacetimes. Using this solution as a bulk spacetime, we can obtain a more smooth 
bouncing universe on the brane. This simple construction works only for k = —1. 

We consider a timelike dust fiuid, of which energy-momentum tensor is given by 

Tf„y = pu^Uy, (3.39) 

where and p are the five-velocity of a fiuid element and energy density, respectively. For 
this matter field in Einstein-Gauss-Bonnet gravity, there is the following exact solution for 
^ [27]: 

ds' = -d? + dx' + r{xf [dx^ + fk{xf{de^ + sin^ e#')] , (3.40) 

P^^^^. (3.42) 

d d 

= ^. (3.43) 

^x^' dt ^ ' 

The solution with = — 1 and A < represents a wormhole supported by a dust fluid with 
negative energy density, of which throat radius is given by 




The quasilocal mass [381 ES] of this solution with A; = — 1 is obtained as 



3(3 + 4aA)\/r^^^,,/,/ 

"5 



It is seen that the energy density of dust is negative, and hence, the weak energy condition 
is violated. Also, the quasilocal mass is everywhere negative. 



This solution can be attached to the exterior vacuum solution f l2.6l) in the GR branch with 
k = —1, A < 0, and a certain negative mass parameter p aX r = r^, where rg is the surface 
radius of the dust region '"s > and h{rs) > are required then. In the resulting 
spacetime, there is no singularity at all. The Penrose diagram of one regularized bulk 
spacetime is drawn in Fig. 3(b) as an example. In summary, the branch singularity for 
k = —1 can be regularized by this matter field. Unfortunately, this simple construction 
does not work in other cases. 
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4 Weakness of the branch singularity 



In this section, we show that the branch singularity is a weak singularity. The physical 
consequence of this property is that a finite body may reach there safely. 

There are several definitions of the strength of a singularity. (See [10|BT|.) We first present 
the definition by Tipler Let 7 : [Aq, Ag) — )• M be an affinely parametrized causal 

geodesic which approaches a singularity as A — )■ A^, where A is an affine parameter. Define 
Jai(7) for Ai G [Ao, Ag) to be a set of maps : [Aq, Ag) — > TM (TM means the tangent 
bundle and / = 1, 2, 3, 4 (/ = 1, 2, 3) for timelike (null) 7) satisfying the following four: 



where A;'" is the tangent of 7. Equation (14. 3 p is called the Jacobi equation (or geodesic 
deviation equation). Along a timelike geodesic, four independent Jacobi fields define a 
volume element V{X) along 7 by the exterior product. Along a null geodesic, three such 
fields define an area element which we also denote V{X). A singularity is called Tipler 
strong if 



is satisfied for all Ai G [Aq, As) and all four (three) linearly independent Jacobi fields Z G 
J\i (7) |12] ■ The singularity is called Tipler weak if it is not Tipler strong. This definition of 
the Tipler strong singularity intuitively says that any object that hits a strong singularity 
is crushed to zero volume (area). 

The above definition ignores the case where V^(A) blows up in the approach to the singu- 
larity. Also, V^(A) may remain finite overall when some of the elements of J\^{'y) blow up 
but some others converge to zero. These possibilities were pointed out by Nolan |33] and 
Ori [H]. In order to include such situations, Ori defined deformationally strong singular- 
ity [H]. A singularity is called deformationally strong if it is either (i) Tipler strong, or (ii) if 
there exists an element of Jai(7) that has infinite norm for A — t- A~ for all Ai G [Aq, Ag) [M]. 
A singularity is called deformationally weak if it is not deformationally strong. 

A singularity is Tipler weak if it is deformationally weak. Here we show that the branch 
singularity is deformationally weak along radial causal geodesies. The first task for this 
purpose is to clarify the asymptotic behavior of causal geodesies near the singularity. 



Zf,)(A) GT,(,)M, 

zU>^i) =0, 




(4.1) 
(4.2) 
(4.3) 
(4.4) 



lim miV{X) = 



(4.5) 
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4.1 Asymptotic behavior of geodesies 



For the metric fl2.6p . the Lagrangian for a geodesic 7 is written as 



(4.6) 
(4.7) 



where a dot denotes the derivative with respect to the affine parameter A. Here k'^ := x'^ 
is a tangent vector of 7 satisfying 



(4i 



where e is given as and —1 for a null and timelike geodesic, respectively. The metric fl2.6p 
does not depend on t and 0, so that from the Lagrange equation we find two independent 
conserved quantities along a geodesic: 



E :-- 



dL 



h{r)i, 



dL 99 9 
$ := ^ = r^f^sin^i 
90 



(4.9) 
(4.10) 



Hereafter we consider only radial geodesies and set the angler coordinates ^5 and all 
constants, which gives $ = 0. Then, the Lagrange equation for r is given as 



r = -h', 
2 ' 

of which the first integral is given from Eq. (14. 8 P as 

f2 = + eh{r). 



(4.11) 



(4.12) 



Now we obtain the asymptotic behavior of the geodesic close to the singularity. First we 
consider the case of timelike 7 (e = — 1). Using the asymptotic expansion (12.130 . Eq. (14.120 
gives the asymptotic behavior of r(A) near r = as 

r(A)^fl(A-As)+f2(A-As)^ (4.13) 

rl.^E^^.[k-l^y (4.14) 

f, :=^f, (4.15) 

A = As is the affine time corresponding to r = 0. 
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On the other hand, using the asymptotic expansion fl2.14p . the asymptotic behavior of r(A) 
near the branch singularity r = rt is given as 

r(A) ^rb + ri(A-As), (4.16) 



r{:=E' ^t\k^^y (4.17) 
Here A = As is the affine time corresponding to r = rb. 

In the case of null geodesies (e = 0), Eq. f l4.12p gets simphfied and is exactly solved to give 

r(A) =rb±^(A-As). (4.18) 

We set ro = rb and tq = for geodesies terminating in or emanating from the branch 
singularity and the central singularity, respectively. 

4.2 Strength of the branch singularity 

We analyze the strength of the branch singularity. We adopt the similar method in four 
dimensions [13]. (See [45j for the higher-dimensional analysis.) The tangent vector to 7 is 
given by 

+ f(A)— , (4.19) 



h{r) dt dr ' 

where r(A) satisfies Eq. fl4.12p . A set of Jacobi fields Z'^j-^ along 7 satisfies the Jacobi 
equation f l4.3p . Because the Jacobi equation is a linear equation for Z^^^, a basis for the 
Jacobi fields can be found by obtaining all independent Jacobi fields in the radial two-space 
and in the tangential three-space. 

In the tangent three-space, there are three independent parallel propagated (namely k'^'VuV(i) 
0) unit spacelike vectors: 

fi d 1 d 

%)^ = -^' (4.20) 
n^^ = I A (4 22) 
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Their dual one-forms e^I^dx'^ (satisfying e^prj^j-^ = 1) are given by 

e^^^dx'' = rdx, (4.23) 
e(f)dx^ = r/,(x)rf^, (4.24) 
e^^Ux^ = rfkix) sin ed(j). (4.25) 

Hereafter, we discuss the cases where 7 is timehke and null, separately. 
4.2.1 Along radial timelike geodesies 

We first consider the case where 7 is timelike. The dual one-form ef^'^dx^ of k^[d/dx^) 
(satisfying e^^^k'^ = —1) is given by 

e^^^dx'' = -hidt + h'^rdr. (4.26) 

A spacelike unit vector in the radial two-space orthogonal to k^^{d/dx^) is given by 



Q r ■ 



r d d 

which is parallel propagated along 7. Its dual one-form e^/^^dx^ is given by 



e 



We consider the following Jacobi fields along 7: 

d , d 



Wdx^" = -rdt + idr. (4.28) 



^inl^, ■= hi)W<i)7^, = 1'2'3,4), (4.29) 



' '^^)^x^' 

of which the norm is given by l^jy The dual one- forms e\!^ of Z(^j) (/ = 1, 2, 3, 4) are given 

by e|/'* := l(i)e^i^\ If 7 is timelike, they define a volume four- form by exterior product and 
its norm is given by 

nA) = |/(i)||/(2)||/(3)|/(4)|. (4.30) 
The behavior of /(/)(A) is determined by the Jacobi equation (14.31) . 

The Jacobi equation for / = 2, 3, 4 gives 

= - ^/(,), (4.31) 
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where we used Eq. fl4.12p and the fact that r^^-^ is parallel propagated. Using Eq. (14. lip . 
Eq. fl4.3ip can be also written as 



This equation is solved to give 



= - '-hiy (4.32) 



ki)W=rW I Z7TV2^ (4-33) 



. r(A) 

which satisfies /(/)(Ai) = 0. Since r is finite for A = [Ao,As), /(/)(A) is nonzero finite for 
A ^ A-. 

For 1 = 1, the Jacobi equation gives 

=Z^(^-^ -\- j.-i.j.Z^-^^-jk^ -\- ^^^Z^^^k^ 

=/.-v(/;i)-ie/i'7(i)), (4.34) 
=Z[i) + R^ urZl)k'k^ + R"- trtZl,)k'k' 
=E(^/(i)-ie/i'7(i)^, (4.35) 

where we used the fact that rj^^-^ is parallel propagated. Hence, both components give the 
same equation: 

=/(!)- ^/(i), (4.36) 

where we used the derivative of Eq. (14. lip at the last equality. Using Eq. (I4.16p . we obtain 
the asymptotic solution around the branch singularity as 



h)W ^ h)o - + o"A(A - (4.37) 

4ar/ V o 

where /(i)o is a nonzero constant and hence /(i)(A) is nonzero finite for A — t- A^. Since all 
Z(7)(A) are nonzero finite for A — )■ A", the branch singularity is deformationally weak along 
radial timelike geodesies. 



4.2.2 Along radial null geodesies 

Next we consider the case where 7 is null. In this case, the Jacobi fields only in the tangent 
space are relevant in order to define the strength of the singularity. The strength is defined 
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by the norm of the area three-form constructed by e*^ {I = 2, 3, 4): 

= 1^(2)1 1^(3)|/{4)|. (4.38) 

From Eq. fl4.3ip . the Jacobi equation for / = 2, 3, 4 is integrated to give 

^(/)(A)=/(/)i(A-Ai), (4.39) 

where is a nonzero constant. This satisfies /(7)(Ai) = and remains finite for A — ?■ Aj7. 
Hence, the branch singularity is deformationally weak also along radial null geodesies. 



5 Summary and discussions 



In this paper, we have presented a novel scenario for the bouncing universe in the five- 
dimensional braneworld in the framework of Einstein-Gauss-Bonnet gravity. In this sce- 
nario, the branch singularity located at the finite physical radius in the bulk, which appears 
for the negative mass parameter, plays an essential role. We have shown that the branch 
singularity can be considered as a massive thin shell and a three-brane may pass through 
it. As a result, the bouncing universe is realized on the brane. 

The bulk spacetime is extended beyond the branch singularity in the C° sense and then the 
branch singularity is identified as a massive thin shell. From the bulk point of view, the 
three-brane collides with another shell of branch singularity and continues to evolve into 
the extended spacetime. From the brane point of view, the moment of the collision is the 
moment of the bounce. This is not the conventional bounce because the derivative of the 
scale factor is discontinuous at the bounce moment. Therefore, there appears a shell-type 
instantaneous singularity on the brane but the curvature invariants never blow up before 
or after then. This claim is strongly supported by the fact that the branch singularity is 
radially deformationally weak, which implies that the singularity is harmless for a finite 
body moving radially. The present result opens a completely new possibility to achieve the 
bouncing brane universe as an effect of the higher-curvature terms. Our scenario is not 
sensitive about the equation of state for the matter on the brane and does work even with 
ordinary matter. 

Here we should make a brief comment on the anisotropy of the universe. Although the 
assumption of the isotropy of the universe is justified at the recombination era by the 
distribution of the cosmic microwave background in the WMAP-Planck data, its validity 
is a nontrivial problem in the very early universe before that. In the standard FRW 
cosmology, it is expected by the BKL (Belinsky-Khalatnikov-Lifshitz) conjecture that such 
an assumption is no more valid and the behavior of the generic spacetime close to the 
initial singularity is quite anisotropic [16]. In contrast, this claim is not necessarily true in 
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the Randall-Sundrum braneworld. It has been shown within the Bianchi class anisotropic 
cosmological models that the brane universe is generically isotropic close to the initial 
singularity. (See [U] for review.) However, the similar analysis has not been done in the 
Gauss-Bonnet braneworld so far. This is one of the important future subjects to clarify the 
validity of our bouncing scenario. 

Our result also gives a suggestion on what would happen if the brane approaches the inner 
Cauchy horizon in the generalized Reissner-Nordstrom-AdS bulk spacetime with a U(l) 
gauge field in the Randall-Sundrum braneworld Pl [12] • The rigorous analysis in the four- 
dimensional asymptotically flat case suggests that the mass inflation instability transforms 
the inner horizon into a weak singularity. More precisely, the resulting singularity appears 
at a finite physical radius where the metric and its inverse are finite [11]. Therefore, the 
extension is possible beyond this singularity. Then a natural question is whether the brane 
can reach and pass through it or not. Of course, the appearance of the singularity means 
that the theory is no more reliable around there. However, the result in the present paper 
suggests that the three-brane could safely pass through the inner horizon in the bulk. 

Since the existence of the branch singularity stems from the quadratic nature of the theory, 
such a singularity is characteristic and must be quite generic in higher-curvature theories. 
In Lovelock higher-curvature gravity [IB], which contains general relativity and Einstein- 
Gauss-Bonnet gravity as special cases, such a singularity appears in this class of vacuum 
solutions rather generically [19] . Interestingly, the central singularity is totally absent and 
the branch singularity is generic independent of the mass paraeter /i if there is a U(l) gauge 
field in the bulk spacetime in Einstein-Gauss-Bonnet gravity [50| |5T] . Undoubtedly, the 
effect of such singularities in cosmology is an interesting problem and should be investigated 
further. 
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A Geometric quantities 



For the metric fl2.6p . the nonzero components of the Christoffel symbol are given by 

T\.-\hh', rs.^|, r..^-|, (A.1) 
= = = ^, (a.2) 



^\. = -rh. r\e = T\,= ^, (A.3) 
r\e = ~rhfl T\, = -hdJ,, (A.4) 

r%, = ^, r,, = -rhf,smH, (A.5) 
Sill u 

T''^^ = -fk{dMsm'9, T\^ = -smecose, (A.6) 

where a prime denotes the derivative with respect to r. The nonzero components of the 
Riemann tensor are given by 

(A.7) 

-Kh'f.sin'e, (A.8) 

^rh'fi sin^ e, (A.9) 

(A.IO) 

(A.ll) 
(A.12) 
(A.13) 
(A. 14) 
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2 ■ 


JDt 

ete — 
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= --rh' 
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JDT 
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-rh'f. 






T3X 

ttx 


- R^ 

— ^ tte - 


td4> 

~ tt4> ~ ' 


hh' 






rrx 


- R^ 


— R"^ — 

^ rr4> 


h' 

2rh' 






TDX 

^ exo 


= fkik- 






- h) sin^ 




^ xxe 


T34> 

~ ^ xx<t> 


= -k + h, 








^ 696 


= fl{k- 


h) sin^ e, 


696 — ' 




h). 



B Asymptotic behavior for a — )► (X) 

In this Appendix, we discuss the asymptotic behavior of the brane universe for a — > oo. In 
the standard Friedmann cosmology f l2.29p . we obtain 

lim V;tand(a) = ^A4. (B.l) 

a— >oo o 

Hence, the recoUapse occurs for A4 < 0. 
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In the Randall- Sundrum braneworld f l2.28p . we obtain 



K, 1 

lim Vcnia) ~ -fa^ + -A 



18a37 a2/ 36a67 Sa^ 



Hence the recollapse condition is given by 



< -6A. 



(B.2) 



(B.3) 



If a precise relation between A and a is satisfied as A = — Kgcr^/G, the condition becomes 
more complicated. 

In the Gauss-Bonnet braneworld with a > 0, the recollapse condition is modified as 
1 



lim VGB(+)(a 



4 \ 3/2 N 1/3 



+ (l + iaA) { (l + laAy\aa'4 + aaV.^l + + 5"a) 

(B.4) 



\ 3/2 



< 0. 



On the other hand, in the case of a < 0, we obtain 



lim VGB(-)(a) 



acT^Kg — acr^Kg^/ 1 + 



4 



+ (l + iaA){-(l + iaA) 



< 0. 



' ( 4 \3/2>|l/3 

qO-H } 

(B.5) 



The above limit is not real if inside the square root is negative. This condition is given by 

2 



am 



4 x3/2 

l + -aAj . 



(B.6) 



If the above inequality is satisfied, there is an upper bound a = Osing in the physical domain 
of CI. CI — '^sing IS db curvature singularity where a and d are finite but d blows up. In the 
case where 



4 



a4\ 3 J 



(B.7) 



is satisfied, the condition (IB. 5 II provides the recollapse condition. 



30 



References 



[1] R. Penrose, Phys. Rev. Lett. 14 , 57 (1965), 

S.W. Hawking, Proc. Roy. Soc. London A300, 187 (1967), 

S.W. Hawking and R. Penrose, Proc. Roy. Soc. London A314, 529 (1970). 

[2] S.W. Hawking and G.F.R. Ellis, "The Large Scale Structure of Space-time" (Cambridge 
University Press, Cambridge, England, 1973); 

F.J. Tipler, C.J.S. Clarke, and G.F.R. Ellis, in General Relativity and Gravitation: 
One Hundred Years after the Birth of Albert Einstein, A. Held, ed. (Plenum Press, 
New York, 1980); 

J.M.M. Senovilla, Gen. Rel. Grav. 30, 701 (1998). 
[3] M. Gasperini and G. Veneziano, Astropart. Phys. 1, 317 (1993). 

[4] M. Novello and S.E.P. Bergliaffa, Phys. Rept. 463, 127 (2008). 

[5] L. RandaU and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999); 
L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 

[6] R. Maartens and K. Koyama, Living Rev. Rel. 13, 5 (2010); 

P. Brax, C. van de Bruck, and A.-C. Davis, Rept. Prog. Phys. 67, 2183 (2004). 

[7] P. Kraus, JHEP 9912, Oil (1999); 
D. Ida, JHEP 0009, 014 (2000). 

[8] P. Hofava and E. Witten, Nucl. Phys. B460, 506 (1996); 
P. Hofava and E. Witten, Nucl. Phys. B475, 94 (1996). 

[9] C. Csaki, J. Erlich, and C. Grojean, Nucl. Phys. B604, 312 (2001); 
C. Barcelo and M. Visser, Phys. Lett. B482, 183 (2000); 
S. Mukherji and M. Peloso, and Phys. Lett. B547, 297 (2002); 
P. Kanti and K. Tamvakis, Phys. Rev. D 68, 024014 (2003). 

[10] E. Poisson and W. Israel, Phys. Rev. D41, 1796 (1990). 

[11] M. Dafermos, Commun. Pure Appl. Math. 58, 0445 (2005); 

M. Dafermos and I. Rodnianski, Invent. Math. 162, 381 (2005). 

[12] J.L. Hovdebo and R.C. Myers, JCAP 0311, 012 (2003). 

[13] N. Okuyama and K.-i. Maeda, Phys. Rev. D 67, 104012 (2003). 

[14] N. Okuyama and K.-i. Maeda, Phys. Rev. D 70, 064030 (2004). 

[15] Y. Shtanov and V. Sahni, Phys. Lett. B557, 1 (2003). 

[16] S. Mukherji and S. Pal, Mod. Phys. Lett. A25, 35 (2010). 



31 



[17] D. J. Gross and E. Witten, Nucl. Phys. B277, 1 (1986); 

D. J. Gross and J. H. Sloan, Nucl. Phys. B291, 41 (1987); 

R. R. Metsaev and A. A. Tseytlin, Phys. Lett. B 191, 354 (1987); 

B. Zwiebach, Phys. Lett. B 156, 315 (1985); 

R. R. Metsaev and A. A. Tseytlin, Nucl. Phys. B293, 385 (1987). 
[18] S.C. Davis, Phys. Rev. D 67, 024030 (2003); 

[19] E. Gravanis and S. Willison, Phys. Lett. B562, 118 (2003); 

C. Garraffo, G. Giribet, E. Gravanis, and S. Wilhson, J. Math. Phys. 49, 042502 
(2008). 

[20] N. Deruelle and C. Germani, Nuovo Cim. 118B, 977 (2003); 

C. Barcelo, C. Germani and C.F. Sopuerta, Phys. Rev. D 68, 104007 (2003). 

[21] C. Charmousis and J.-F. Dufaux, Class. Quant. Grav. 19, 4671 (2002). 

[22] J.E. Lidsey and N.J. Nunes, Phys. Rev. D 67, 103510 (2003); 

G. Kofinas, R. Maartens, and E. Papantonopoulos, JHEP 0310, 066 (2003); 
B.C. Paul and M. Sami, Phys.Rev. D 70, 027301 (2004); 

X.-He Meng and P. Wang, Class. Quant. Grav. 21, 2527 (2004); 
M. Sami and V. Sahni, Phys. Rev. D 70, 083513, (2004); 
A. Padilla, Class. Quant. Grav. 22, 681 (2005); 
T. Kobayashi, Gen. Rel. Grav. 37, 1869 (2005); 

R.A. Brown, R. Maartens, E. Papantonopoulos, and V. Zamarias, JCAP 0511, 008 
(2005); 

R.A. Brown, Gen. Rel. Grav. 39, 477 (2007); 
K. Konya, Class. Quant. Grav. 24, 2761 (2007); 
K. Konya, Phys. Rev. D 75, 104003 (2007); 

M. Heydari-fard and H.R. Sepangi, Phys. Rev. D 75, 064010 (2007); 

P.S. Apostolopoulos, N. Brouzakis, N. Tetradis, and E. Tzavara Phys. Rev. D 76, 

084029 (2007); 

H. Maeda, V. Sahni, and Y. Shtanov, Phys.Rev.D 76, 104028 (2007), Erratum-ibid. 
D 80, 089902 (2009); 

S.-F. Wu, A. Chatrabhuti, G.-H. Yang, and P.-M. Zhang, Phys. Lett. B659, 45 (2008); 
S. del Campo, R. Herrera, and J. Saavedra, Mod. Phys. Lett. A23, 1187 (2008); 

E. N. Saridakis, Phys. Lett. B661, 335 (2008); 

M. Bouhmadi-Lopez and P.V. Moniz Phys. Rev. D 78, 084019 (2008); 

K. Nozari and B. Fazlpour, JCAP 0806, 032 (2008); 

K. Nozari and N. Rashidi, Int. J. Theor. Phys. 48, 2800 (2009); 

K. Nozari and T. Azizi, Phys. Scripta 83, 015001 (2011); 

R. Herrera and N. Videla, Eur. Phys. J. C67, 499 (2010). 

[23] R.A. Brown, ^^Brane world cosmology with Gauss-Bonnet and induced gravity terms" 
(PhD thesis. University of Portsmouth, 2006), [ir^/0701083 



32 



[24] K.-i. Maeda and T. Torii, Phys.Rev. D 69, 024002 (2004); 

T. Kobayashi, T. Shiromizu, and N. Deruelle, Phys. Rev. D 74, 104031 (2006). 

[25] S. Mukohyama, Phys. Rev. D 65, 084036 (2002); 

N. Deruelle and M. Sasaki, Prog. Theor. Phys. 110, 441(2003); 
G. Calcagni, Phys. Rev. D 69, 103508 (2004); 

M. Minamitsuji and M. Sasaki, Prog. Theor. Phys. 112, 451 (2004); 

S. Tsujikawa, M. Sami, and R. Maartens, Phys. Rev. D 70, 063525 (2004); 

S.C. Davis, Phys. Rev. D 72, 024026 (2005); 

R.-G. Gai, H.-S. Zhang, and A. Wang, Gommun. Theor. Phys. 44, 948 (2005); 
G. Panotopoulos, Nucl. Phys. B745, 49 (2006); 

J.-H. He, B. Wang, and E. Papantonopoulos, Phys. Lett. B654, 133 (2007); 

A. Herrera-Aguilar, D. Malagon-Morejon, R.R. Mora-Luna, and I. Quiros, 

Glass.Quant.Grav. 29, 035012 (2012). 

[26] S. Nojiri, S.D. Odintsov, JHEP 0007, 049 (2000); 

K. Aoyanagi and K.-i. Maeda, Phys. Rev. D 70, 123506 (2004); 
S. Nojiri and S.D. Odintsov, Gen. Rel. Grav. 37, 1419 (2005). 

[27] J.-F. Dufaux, J.E. Lidsey R. Maartens, and M. Sami, Phys. Rev. D 70, 083525 (2004); 
M. Sami, N. Savchenko, and A. Toporensky, Phys. Rev. D 70, 123528 (2004); 
D. Konikowska, M. Olechowski, and M.G. Schmidt, Phys. Rev. D 73, 105018 (2006); 
T. Kobayashi and M. Minamitsuji, JGAP 0612, 008 (2006). 

[28] LP. Neupane, JHEP 0009, 040 (2000); 

Y.M. Gho and LP. Neupane, and P.S. Wesson, Nucl. Phys. B621, 388 (2002); 
Y.M. Gho and LP. Neupane, Int. J. Mod. Phys. A18, 2703 (2003); 
G. Gharmousis, S.G. Davis, and J.-F. Dufaux, JHEP 0312, 029 (2003); 
Ph. Brax, N. Ghatillon, and D.A. Steer, Phys. Lett. B608, 130 (2005). 

[29] N.E. Mavromatos and J. Rizos, Phys. Rev. D 62, 124004 (2000); 
N.E. Mavromatos and J. Rizos, Int. J. Mod. Phys. A18, 57 (2003); 
P. Binetruy, G. Gharmousis, S.G. Davis, and J.-F. Dufaux, Phys. Lett. B544, 183 
(2002); 

J.P. Gregory and A. Padilla, Glass.Quant.Grav. 20, 4221 (2003). 

[30] T. Glifton, P.G. Ferreira, A. Padilla, and G. Skordis, Phys. Rept. 513, 1 (2012). 

[31] D. G. Boulware, and S. Deser, Phys. Rev. Lett. 55, 2656 (1985). 

[32] J. T. Wheeler, Nucl. Phys. B268, 737 (1986); 

D. Lorenz-Petzold, Mod. Phys. Lett. A3, 827 (1988); 

R.-G. Gai, Phys. Rev. D 65, 084014 (2002); 

R.-G. Gai and Qi Guo, Phys. Rev. D 69, 104025 (2004). 

[33] R.M. Wald, General Relativity, (University of Ghicago Press, 1984). 



33 



C. Charmousis and A. Padilla, JHEP 0812, 038 (2008). 
T. Torii and H. Maeda, Phys. Rev. D71, 124002 (2005). 

E. Gravanis and S. Willison, Phys. Rev. D 75, 084025 (2007). 
H. Maeda and M. Nozawa, Phys. Rev. D 78, 024005 (2008). 
H. Maeda, Phys. Rev. D 73, 104004 (2006). 

H. Maeda and M. Nozawa, Phys. Rev. D77, 064031 (2008). 

P.S. Joshi, Global Aspects in Gravitation and Cosmology (Oxford University Press, 
New York, 1993). 

C.J.S. Clarke, The analysis of Space- Time Singularities (Cambridge University Press, 
Cambridge, 1993). 

F. J. Tipler, Phys. Lett. A 64, 8 (1977). 

B.C. Nolan, Phys. Rev. D 60, 024014 (1999); 
B.C. Nolan, Phys. Rev. D 62, 044015 (2000); 

A. Ori, Phys. Rev. D 61, 064016 (2000). 

S. Abdolrahimi and A. A. Shoom, Phys. Rev. D 81, 024035 (2010). 

J. Wainwright and G.F.R. Ellis, "Dynamical Systems in Gosmology" (Cambridge Uni- 
versity Press, Cambridge, England, 1997). 



A.A. Coley, e-Print: |astro- ph/0504226. 
D. Lovelock, J. Math. Phys. 12, 498 (1971). 

H. Maeda, S. Wilhson, and S. Ray, Class. Quant. Grav. 28, 165005 (2011) 

D.L. Wihshire, Phys. Lett. B169, 36 (1986); 
M. Cvetic, S. Nojiri, S.D. Odintsov, Nucl. Phys. B628, 295 (2002). 

[51] T. Torii and H. Maeda, Phys. Rev. D72, 064007 (2005). 



34 



